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Abstract 

An estimate has been made of the masses of heavy hadrons in nonrelativistic quark 
model, which includes spin and flavor-dependent hyperfine splitting for two quarks. 
The effect of variation of the wavefunction value at origin i.e. V(0)| 2 and the strong 
coupling constant a s , with flavor, has also been included in calculating the mass values. 
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1 Introduction 



A lot of data are available on the masses of mesons but the masses of most baryons have 
not been measured yet. Recently predictions about the heavy baryon mass spectrum have 
become a subject of increasing interest due to the current experimental activity of several 
groups at CERN, Fermilab and the Cornell Electron Storage Ring (CESR), aimed at the 
discovery of the baryons so far absent from the baryon summary table [1]. The copious 
production of heavy quarks at LEP, Fermilab Tevatron, CERN LHC and B factories, open 
for study the rich spectroscopy of heavy hadrons. So, a plausible theoretical prediction for 
the baryon mass spectrum becomes a guide for experimentalists. 

Several models [2-10] have been used to evaluate the heavy baryon mass spectrum. The 
nonrelativistic quark model (NRQM), simple and economic one, has been able to explain 
very nicely the mass spectrum of light baryons and mesons. Many workers [3,4] have studied 
the masses of heavy hadrons in NRQM. Long time ago, Singh and Khanna [3] , ignoring the 
effect of variation of strong coupling constant a s and IVKO)! 2 (the wave function value at 
origin) with flavor, had estimated the masses of heavy hadrons using the NRQM with the 
inclusion of spin and flavor-dependent hyperfine interaction between two quarks and between 
a quark and an antiquark. They assume IV^O)! 2 scale for the heavy baryons to be the same 
as that of hyperons. However, since IV^O)! 2 is a dimensional quantity it may be incorrect to 
ignore its variation with flavor. Evidence to corroborate this has been found in the quark 
model [11-13] as well as in the lattice calculation [14]. Moreover, a s is also scale dependent, 
so it will be worthwhile to include its variation with the mass scale. 

In this paper, we calculate the masses of heavy hadrons in NRQM with the inclusion of 
variation of IVKO)! 2 an d ot s with flavor. 

2 Mass Formulae 

The mass of a heavy hadron is assumed to arise from the constituent quark masses plus 
the two-body hyperfine interaction energies for a meson and a baryon. Without taking 
into consideration the variation of IVKO)! 2 and a s with flavor, hadron masses, thus, can be 
expressed as [3], 

s* s * 

M moson = a m +mi+ rrij + bm Q '" J , (1) 

niiinj 
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M baryon = « b +E^ + ^rE5S/ P) 

i i<j J 

where a m and on, are the parameters with the dimension of mass and m, (mj ) are the masses 
of respective quarks (antiquarks). m is a scale factor which is taken to be the mass of the 
lightest quark i.e. m = m u . 'b', having the dimensions of mass, is a parameter which 
includes the long distance effects i.e. it takes care of the value of the wavefunction at origin, 
IV^O)! 2 , along with the value of strong coupling constant a s and is usually taken to be the 
same for all hadrons, irrespective of the quark flavor. 

Now, if the dependence of 1^(0) | 2 and a s on flavor is taken into account, eqs. (1) and 
(2) modify to 

s* s * 

M mcson = a m + mi + mj + 6i,-m '" J , (3) 

nurrij 

2 

M baryon = a b + J2 m < + ^-J^bij-^-. (4) 

i %<j J 

Here, the scale factor, b^, instead of having a constant value for all flavors, is a variable and 
thus, has different values for different quark pairs. Hence, b^ is the parameter which takes 
into consideration the flavor dependence of \tp(0)\ 2 and a s , so that one has, 

hj oc \ipij(0)\ 2 a s (n). (5) 

The value of bij is taken to be same in the light quark sector i.e. for any of the u, d and s 
quark pairs and it varies as one goes to c-sector and then to b-sector. 



3 Estimation of |^(0)| 2 and a s (/i) 
3-1 |^(0)| 2 

The flavor dependence reflected in the scale factor bij corresponding to the spatial matrix 
element, is due to the long distance QCD effects. Evaluation of IV^O)! 2 is as yet uncertain 
for baryons and more complicated, because unlike the mesons, these are three-body systems. 
Infact, the heavy baryons may provide a good way of testing the flavor dependence of the 
confinement forces. The absence of an exact dynamical theory of low energy interactions 
between quarks limits our evaluation of |^(0)| 2 from first principles. However, a naive 
estimate for the scale parameter may be obtained using the hyperfine splitting (hfs) [2], 

1 f\ IT f\' 

AE hfs = ^-\m\^,s h (6) 
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The experimental hyperfine splittings thus, may provide a reliable measure of the wave- 
function at origin of the ground state baryons the value of which is needed in the lifetime 
estimates. 

Many workers have tried to estimate the value of |?/>(0)| 2 for different quark flavors 
using different techniques [13-18]. H. Y. Cheng [15], using formula (5) has found the values 
respectively, for c and b sectors to be 

|<M0)| 2 = 1.5xl(T 2 GeV 3 (7) 
and 

|^„(0)| 2 = 2.5 x lCr 2 GeV 3 . (8) 

Korner and Siebert [16], from a fit to hyperfine splitting, have estimated the value of IV^O)! 2 
for c-sector to be nearly 1.0 x 10 -2 GeV 3 . As mentioned in ref. [16], NRQM with a funnel 
type potential [17] and electromagnetic mass differences in the hyperfine splitting formula 
[18] also predict the values similar to that found by them. Values of the scale parameter 
IV^O)! 2 for different quark pairs have also been estimated in lattice approach and have been 
quoted in ref. [14]. 

Uppal and Vcrma [13], using cqn. (5) have made an estimate of the ratio |^[q]|s , which 
comes out to be approximately 2.83. We have also tried to estimate the wavefunction value 
at origin in terms of the ratios, |V'cm| 2 /|V's«| 2 and |V'6«| 2 /|V's«| 2 , respectively, for c and b 
sectors. 

3.2 a a (ii) 

The QCD coupling constant a s (/z) at any rcnormalization scale can be calculated from 
a s (raz) = 0.117 via 

asW ~ l-(ll-fn / )[a a (mz)/27r]/n(mz//i)' U 

and one has, 

a s (m c ) = 0.31 and a s (m b ) = 0.20 (10) 



4 Numerical Estimate for Masses 
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4.1 Heavy Baryons 

Choosing the set of parameters (in MeV) to be a& = 205, m = m u = rrid = 310, m s = 
450, m c = 1620, m b — 4976, b qq < — 640, b qc = 736 and b qb = 1600, with q and g' = u, d, s, 
the masses of heavy baryons are calculated. Following the usual convention that a particle 
symbol represents its mass, the results for c-sector are displayed in Table 1 whereas Table 
2 contains b-baryon masses. The mass values are in good agreement with the experimental 
values, wherever available. For comparison, masses estimated by other techniques have also 
been included in respective Tables 1 and 2. 

Predictions have also been made for the masses of doubly heavy baryons using \ip qc (0) | 2 = 
|V'cc(0)| 2 and \ip q b(0)\ 2 — | V'/ifc (0) | 2 , with h — b or c. Although it may not be an appropriate 
assumption, as no data are available in this sector, we can still make this choice. Moreover, 
when a particle contains two heavy quarks, in the heavy quark mass limit, the corresponding 
hyperfine splitting term (oc l/m,mj) will give negligible contribution to the overall mass 
and hence will be less significant. Predictions made on this assumption along with those 
made by other workers, have been given in Tables 1-3. 

4.2 Heavy Mesons 

Many different sets [8,19-29] of constituent quark masses to evaluate the masses of heavy 
mesons, are used in the literature, most of them obtained from fits to spectroscopic data. 
The values (in MeV) of various parameters used by us are: a m = 82, m = m u = m d = 
310, m s = 400, m c = 1578, m b = 4920, b qq , = 645, b qc = 684 and b qb = 755, with 
q and q' — u, d, s. Mass values thus obtained are displayed in Table 4 and arc in nice 
agreement with the experimental values. 

It is interesting to note that the quark masses which give a best fit to the baryons are a 
little higher (except for the mass of u quark) than those which lead to a best fit to the mesons. 
Because these quark masses are constituent quark masses, there are no theoretical reasons 
why the masses determined from the baryons should coincide exactly with those determined 
from the mesons. If we insist that a single set of quark masses hold for both baryons and 
mesons, and vary these masses, the overall best fit to the hadron data is significantly poorer 
and our predictions have greater errors. 
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5 Discussion 



At present, many of the charm baryon masses are known experimentally, whereas in b-sector, 
only the mass of one b-baryon is known accurately. The charm baryon masses (in McV) 
measured till date are [1] 

A c = 2285, 
S c = 2453 ±1, 
E c = 2468 ± 2, 
Q c = 2704 ± 4, 
E* = 2521 ±4, 

E* = 2644 ±2. (11) 

In 1994, WA89 Collaboration [30] had reported the observation of E' c = 2563 ± 15 McV. 
Recent calculations by many workers [8-10,31-34] consistently predict the mass of E' c to be 
around 2580 MeV. Also, the preliminary CLEO results on the E' c — E c mass difference from 
E t' -> s ^7, -> S°7 decay [35] give us: 

S+'-S+ = 107.8 ± 1.7 ±2.5MeV, 

E Q C '-E a c = 107.0 ± 1.4 ± 2.5MeV. (12) 

With E c = 2468 ± 2 McV, these results lead to 

E' c = 2575 ± 5MeV. (13) 

Our prediction for E' c (2580 MeV) is in nice agreement with this experimental observation as 
well as with the predictions made by other authors. Also the mass difference E' c — E Cl using 
our mass values, comes out to be 107 MeV, which coincides very well with the experimentally 
measured value. 

The f2* as well as double and triple charm baryons, have not yet been observed. Most of 
the theoretical calculations [7-10,31-34,36] predict the mass of 17* to be around 2770 MeV, 
whereas predictions for the doubly and triply charmed baryon masses are less definite. The 
mass value, as predicted by us, for fi* (2766 MeV) is consistent with the above mentioned 
theoretical predictions. 

On the other hand, in the b-sector, only the mass of A& [1] baryon is known well and 
masses of the rest of the b-baryons are yet to be measured. The mass of A& baryon as 
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predicted by us is in nice agreement with the experimental value [1], whereas the rest of 
b-baryon masses predicted by us agree well with the predictions made by ref. [8] , as can be 
seen from Table 2. 

To have an estimate of the wavefunction value at origin, we make use of eqn. (5) i.e. 
bij oc | VYj (0)| 2 c*s(m)> s ° that one has 

|^c«(0)| 2 a s (m c ) 



\tp su (0)\ 2 a s (m s ) 

and 

hu \ipbu(0)\ 2 a s (m b ) 
b S u \4>su(0)\ 2 a s (m s ) 

Using our preferred set of parameters, these ratios turn out to be 



(14) 



(15) 



and 



= 1.15 (16) 



^ = 2.5 (17) 



As the mass dependence of the strong coupling constant (cqn. 9) is known, an estimate of 
the ratios of the wave function values at origin (i.e. I^cmP/I^smI 2 and |V>6«| 2 /IVVu| 2 ) can be 
made using eqs. (15) to (18): 

|Vw|7lV> s «| 2 = 3.3 
and |Vb„| 2 /IVw| 2 = 11.3 (18) 

The ratio for c-sector is little higher than the estimates of ref. [13] and that of lattice 
calculations [14]. It may be due to the fact that the ratio also depends on the strong 
coupling constant for light sector, whose value is not fixed. So different values of a s (m s ) 
will lead to different values of these ratios. Only the experimental data in future will give 
us better insight in this arena. 

Note that, irrespective of introducing the parameters concerning the variation of wave- 
function at origin with flavor, certain independent relations and mass shifts introduced in 
ref. [3] are satisfied here also: 

fl* b + - 2E* b = Q* c + S* - 2S* = fl* + S* - 2S*, (19) 

(E* c - K) = S cc - Sc, (20) 
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- £*) = E bb - E 6) (21) 

(s; - gb) = (s; - gc) = (s* - g) = A-jv = i 

(B* — B) (D* -D) (K* -K) p-n 2 [ ' 

3 S b - A fc B* - B 3 £ c - A c £>*-£> 3 £ - A K* - K 
2A-N p-ir 2A-N p-n 2A-N p-n 1 ; 

The predicted masses may have large errors, so it may be useful to have an idea of the mass 
difference from the lowest lying heavy baryon. In b-sector, one has (in MeV): 





A b = 


180, 


£*- 


A b = 


230, 


— 


A 6 = 


188, 




A 6 = 


302, 


— * 


A 6 = 


343, 




A 6 = 


436, 


n* b - 


-A b = 


= 469 



(24) 



To conclude, by allowing the variation of parameter b with flavor and introducing two more 
parameters, we are able to bring lot of masses of both baryons and mesons closer to the 
observed values. 
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Table 1: Masses (in MeV) of charm baryons 



Particle 


Calculated Mass 


Experimental Value 


Rcf. [3] 


Rcf. [8] 


Ref. [34] 


A c 


2285 


2285 


2267 


2285±1 


- 


s c 


2451 


2455 


2436 


2453±3 


- 


1 — 'C 


2475 


2470 


2504 


2468±3 


- 


1 — 'c 


2582 


2580* 


2602 


2580±20 


2569±6 


«c 


2717 


2704 


2775 


2710±30 


- 


£c* 


2521 


2521 


2491 


2520±20 




"* 


2642 


2644 


2646 


2650±20 




ft* 


2766 




2800 


2770±30 


2767±7 




3714 




3710 


3660±70 


3610±3 


"cc 


3786 




3781 


3740±70 


3735±17 


ftcc 


3882 




3865 


3740±80 


3804±8 


n* 


3914 




3948 


3820±80 


3850±25 



Table 2: Masses (in MeV) of beauty baryons 



Particle 


Calculated Mass 


Experimental Value 


Rcf. [3] 


Ref. [8] 


A fc 


5641 


5641 


5724 


5620±40 


s 6 


5821 


5814* 


5729 


5820±40 




5831 




5783 


5810±40 




5949 




5893 


5950±40 


n b 


6083 




6065 


6060±50 


£6* 


5871 


5870* 


5750 


5850±40 




5991 




5915 


5980±40 


ft* 


6117 




6079 


6090±50 


2 bb 


10434 




10315 


10340±100 




10484 




10321 


10370±100 


ft&6 


10585 




10494 


10370±100 


ft&6 


10619 




10321 


10400±100 



* as calculated from the mass differences. 
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Table 3: Masses (in MeV) of baryons containing two different heavy quarks i.e. c and b. 



Particle 


"6c 


"be 


"* 

"fcc 








Mass 


7076 


7102 


7133 


7227 


7244 


7266 


Rcf. [8] 


6990±90 


7040±90 


7060±90 


7060±90 


7090±90 


7120±90 



Table 4: Masses (in MeV) of heavy mesons. 



Particle 


Calculated Mass 


Experimental Value [1] 


D 


1869 


1869 


D* 


2003 


2007 


D s 


1982 


1968 


D* 


2086 


2112 


B 


5276 


5279±1.8 


B* 


5324 


5325±1.8 


B s 


5374 


5369±2 


B* s 


5411 
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